Critical exponents of the quark-gluon bags model with the critical endpoint 
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The critical indices a' , f3, 7' and S of the Quark Gluon Bags with Surface Tension Model that has 
the critical endpoint are calculated and compared with the exponents of other models. These indices 
are expressed in terms of the most general parameters of the model. Despite the usual expectations 
the found critical indices do not depend on the Fisher exponent r and on the parameter x which 
relates the mean bag surface to its volume. The scaling relations for the obtained critical exponents 
are verified and it is demonstrated that for the standard definition of the index a' the Fisher and the 
Griffiths scaling inequalities are not fulfilled in general case, whereas the Liberman scaling inequality 
is always obeyed. This is not surprising for the phase diagram with the asymmetric properties of 
pure phases, but the present model also provides us with the first and explicit example that the 
specially defined index a' s does not recover the scaling relations as well. Therefore, here we suggest 
the physically motivated definition of the index a' = a' c and demonstrate that such a definition 
recovers the Fisher scaling inequality, while it is shown that the Griffiths inequality should be 
generalized for the phase diagram with the asymmetric properties. The critical exponents of several 
systems that belong to different universality classes are successfully described by the parameters of 
the present model and hence its equation of state can be used for a variety of practical applications. 
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I. INTRODUCTION 

Scaling has been widely accepted as "a pillar of mod- 
ern critical phenomena" p] . The scaling hypothesis used 
in the study of critical phenomena was independently de- 
veloped by the well-known scientists as Widom, Domb, 
Hunter, Kadanoff, Fisher, Patashinskii and Pokrovskii 
(see review [2] for the details). Further it was developed 
within the renormalization group approach [5J 0] . One of 
the most striking predictions of the scaling hypothesis is 
called the scaling laws. For instance, for the ordinary liq- 
uids these scaling laws relate the critical exponents a, ft, 
7 and 5 which describe the behavior of the specific heat 
capacity (C ~ |i|~ Q \ here t — T y T " p < denotes a rela- 
tive deviation of the temperature T from the critical one 
T cep ), density differences of the liquid and gaseous phases 
(pi — P g ~ |*|^), isothermal compressibility (Kt ~ |i|~ 7 ) 
and the shape of the critical isotherm which is given by 
the critical index 5 (for the formal definitions see below) 

Fisher [5] : a' + 2/3 + 7' > 2 , (1) 
Griffiths [6] : a! +$(1 + 8) > 2 , (2) 
Liberman [7] : 7' + - £) >0. (3) 

Similar equalities can be also introduced for magnetic 
systems [8] and for percolating systems The corre- 
sponding exponent inequalities for magnetic systems are 
often called Rushbrooke's [TU], Griffiths' and Widom's 
[llj inequalities, respectively. 

The superscript prime in Eqs. ([l])-([3| is necessary to 
introduce for the systems with the phase diagram for 
which the behavior of specific heat capacity and/or com- 
pressibility on the gaseous side of phase diagram, where 
t < 0, differs from that one on the liquid side, i.e. for 



t > 0. Although in many physical systems the scaling 
laws ([ij— ([3]) are obeyed as equalities, it is customary to 
write them as inequalities since, as was proven by M. 
E. Fisher for liquids [T2], in the most general case 
they can be established as inequalities only, ft was, how- 
ever, found that in such exactly solvable models as the 
Fisher- Felderhof one-dimensional model [13], the statis- 
tical multifragmentation model (SMM) [T3], the quark 
gluon bags with surface tension model (QGBSTM1) with 
the tricritical endpoint [15] and its generalization |16j . 
the Fisher and the Griffiths scaling inequalities ([!]) and 
([2]), in which the index a' is involved, may be broken. 
The corresponding proofs are given in [13 for the Fisher- 
Felderhof model, in [17] for the SMM and in [18] for the 
QGBSTM1 and its generalization [IB]. In all these cases 
the phase diagram is rather asymmetric, i.e. the behavior 
of pure phases on the both sides of the phase equilibrium 
curve are different, and, as a result, the heat capacity of 
gaseous and liquid phases are quite different. To some ex- 
tent this problem was resolved by M. E. Fisher in [T3J by 
introduction of a specially defined index a' s which mea- 
sures a divergency of the heat capacity difference of two 
phases at the critical endpoint (CEP). In this case the 
a' s index gives the maximal value among the a' index of 
the gaseous and the liquid phases. The usage of such an 
index instead of the traditional exponent a' allows one to 
formally recover the Fisher and the Griffiths inequalities 
in all exactly solvable models mentioned above, although 
neither the physical meaning of the index a' s nor its rela- 
tion to the experimental procedure of the heat capacity 
measurement were ever justified. Therefore, we are faced 
to three principal questions: (I) Is it possible to justify 
the a' s definition? (II) What definition should be used 
instead of the a' s index, if it fails to recover the scaling 
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inequalities ([T])-((3|)? (Ill) What should we do with the 
scaling laws in the latter case? 

In order to clarify these questions in the present work 
we calculate the critical exponents of the quark gluon 
bags with surface tension model (QGBSTM2 hereafter) 
which has the CEP [19] . This phenomenological model is 
a novel development of the well known gas of bags model 
[20j which, however, contains entirely new mathemati- 
cal mechanism of the CEP generation. Of course, the 
QGBSTM2 was developed for the CEP modeling with- 
out specifying its universality class and the present work 
is devoted to the finding of the critical exponents of this 
model and to the determination of its universality classes. 
Since our main subject of interest is to describe the end- 
point properties of the deconfinement phase transition 
(PT) of quantum chromodynamics (QCD), we would like 
to pay a special attention to the QCD CEP properties. 
For this purpose we, in accord with the contemporary 
knowledge, suppose that the CEP of the deconfinement 
PT has the properties typical for 3-D Ising model in the 
case of 3 quark flavors degenerated QCD [2TM23] , whereas 
for 2+1 quark flavors we assume that the QCD endpoint 
belongs to the universality class of the 0(4) symmet- 
ric 3-dimensional spin model [2~4Tl28| . In other words we 
would like to determine the QGBSTM2 parameters and 
fix them in order to reproduce the critical exponents of 
respective universality class. 

In contrast to the comparable solvable models [T3TTTS] 
including the QGBSTM1, the QGBSTM2 has entirely 
different structure of isobaric ensemble singularities de- 
scribing the PT, and, as we show below, for some values 
of parameters the Fisher and the Griffiths inequalities (fTl) 
and ([2]) of this model are not fulfilled for both a' and a s 
indices. Therefore, here we introduce a physically moti- 
vated definition of the supremum index a' c which is found 
from the linear combination of the specific heat capacities 
of pure phases taken with the nonsingular weights. On 
the examples of solvable models [TBTHo] we demonstrate 
that such a definition recovers the scaling inequalities 
in those case, when the traditional a' index fails, since 
c^c — a 's- However, below we show that even with an im- 
proved definition of the index a' = a' c for the QGBSTM2 
critical exponents there are two regimes: the traditional 
scaling regime, when the scaling inequalities ([l])-((3]) are 
held as the equalities only, and the generalized scaling 
regime, when the Fisher and the Liberman ([3| in- 
equalities are fulfilled, but the Griffiths inequality (f2j) is 
only obeyed in its generalized form which is suggested 
here. The performed thorough analysis of the scaling 
laws ([l}-((3]) and their generalization onto the case of the 
phase diagram with the asymmetric properties seem to 
be very important nowadays in a view of fast techno- 
logical and computational progress which, respectively, 
allows one to study the substances and models with new 
and unusual thermodynamic properties. 

The work is organized as follows. Section [TT] is devoted 
to a brief discussion of the QGBSTM2 main ingredients. 
The model is analyzed in details in Section ITTT) The QG- 



BSTM2 critical exponents are calculated in that section 
also. The analysis of the scaling relations between the 



found critical exponents is given in Section [TV 
sions and perspectives are discussed in Section 



Conclu- 



QUARK GLUON BAGS WITH SURFACE 
TENSION MODEL 



II. 



The QGBSTM2 |T5l [19] treats the quark-gluon plasma 
(QGP) bags and hadrons as relevant degrees of freedom. 
Similarly to the original statistical bootstrap model [25] 
the attraction between the degrees of freedom in this 
model is accounted via many sorts of the constituents, 
while the repulsion between them is introduced a la Van 
der Waals equation of state [151 HH [2U] . The phase struc- 
ture of the QGBSTM2 is completely defined by the mass- 
volume spectrum that for a given temperature T, bary- 
onic chemical potential /1 is defined as 

F{z,T,fi) = F H {z,T,fi)+u(T,n)I T {Az,T,). (4) 

This spectrum defines the isobaric partition [T5J H3 

1 



Z(z,T,fi) 



z - F(z,T,fi) ' 



(5) 



where z denotes the isobar variable. 

The discrete part of the mass-volume spectrum Fh in 
Q is successfully used as the hadron resonance gas model 
to describe the experimental hadron multiplicities which 
allow one to recover the thermodynamic quantities of 
strongly interacting matter created in the heavy ion col- 
lisions, when this matter reaches the chemical freeze-out 
stage (an incomplete list of related works can be found 
in [30l [31] ) . Here we consider the simplest parameteri- 
zation of Fh, since both the quantum statistics and the 
width of hadron resonances are important for the tem- 
peratures below 50 MeV and the baryonic chemical po- 
tentials larger than 940 MeV [30] [3]]. Thus the spectrum 
Fh is parameterized as follows 

n 

F H (z,T,ix) = ^ 5i e^-»'^(T, mi ). (6) 
i=i 

The continuous part of the spectrum Q is chosen in the 
simplest form (compare it that one used in |16j). which 
can be cast as an integral (Vq ~ 1 fm 3 [3"2"] ) 



J T (Az,E) = 



dv 



(7) 



v 



where the Fisher exponent r > 2 provides the conver- 
gence of the integral §f§ for Az — and £ — 0. Also 
here the notation Az = z — Zm(T,/j,) is introduced. In 

tthe particle density of a hadron of mass rrij , baryonic 
irge bj, eigenvolume Vj and degeneracy gj is denoted 
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as 



(j>(T, rrij 



1 

2^ 



p dpe 



(p 2 



(8) 



The functions u(T, /x) and zjy(T, /x) in Q and ([7| are the 
parameters of the present model which are smooth and 
finite together with all their first and second derivatives 
[T51 [T§] . The z-dependent exponentials in ([6| and Q 
describe the short range repulsion of the Van der Waals 
type [13 US] . To parameterize the surface of a QGP bag 
in the continuous part of the spectrum ^ the parameter 
k is introduced. Usually the constant >c is defined by the 
dimension d as x = , but in what follows it is treated 
as a free parameter with the range of values < x < 1. 

A few words should be added here about the hadronic 
surface tension. In principle, it can be included into the 
QGBSTM2 discrete spectrum The first and inter- 
esting results about the surface tension of hadrons which 
fit well into the QGBSTM2 framework can be found in 
[3"T] . They clearly demonstrate that the hadronic surface 
tension is rather small, although it changes the sign at 
the temperature about 150 MeV. However, in the present 
work we do not consider this element and set the hadronic 
surface tension to zero since its inclusion does not affect 
the expressions for the critical indices. 

The new element of principal importance of the present 
model is the parameterization of the surface tension coef- 
ficient TE(T,n) (E(T,/i) in (7) denotes the reduced sur- 
face tension coefficient) which in the vicinity of the phase 
equilibrium curve T = T c (/i) is defined as 



T - T c (p cep ) + ^-(fj, ce p - (J,) 



dfi 



T 



(9) 



with the following values of constants > 1, ^ ± > 0. 
Here <jq is chosen to be a positive constant, but the ob- 
tained results hold, if ctq > is a smooth function of T 
and [i. As shown below it is also of crucial importance 
that the parameters C and have different values be- 
low and above the phase coexistence curve T = T c (/i) 
which exists for /i > p cep . It can be shown [19] that the 
necessary condition for the deconfinement PT existence 
with the CEP is that the QGBSTM2 surface tension co- 
efficient changes the sign exactly at the phase equilibrium 
curve. In other words, the solution Ts(/i) of the equation 
E(T, /i) = should coincide with the PT curve T c (/j,) for 
|U > fi cep , i.e. T S ((J,) = T c {fi) for p > p cep (see Fig. 1 for 
details). The important physical consequence of the nil 
surface tension curve and the deconfinement PT curve 
matching is that for the corresponding thermodynamic 
parameters the volume (and, hence, the mass) distribu- 
tion of QGP bags becomes the power-like instead of ex- 
ponential. Typically, in the liquid-gas PT the power-like 
volume distribution of droplets corresponds to the CEP 



[3"3"] and it leads to the formation of fractals. Such a 
power law may, in principle, naturally explain the ap- 
pearance of the non-Boltzmann fluctuations in the high 
energy collisions experiments .34, 35 without appealing 
to the Tsallis statistics [55] , 

The important mathematical consequence of such a 
matching is that the discontinuity of the partial \x and 
T derivatives of the reduced surface tension coefficient 
across the line T — Ts(/i) provides the 1st order de- 
confinement PT existence [TH], and hence in ([9| one has 
( + 7^ £~ and £ + ^ £~ in general. The quantities in- 
troduced in (|9| have the superscript +(— ), if they are 
taken for T above (below) the curve Ts(m) m the whole 
fx — T plane. For < fi < \i cep the nil line of the surface 
tension coefficient is located in such a way that the de- 
confinement PT degenerates into a cross-over since in this 
region E(/z) < (see [TH] and Fig. 1 for more details). 

Also note that the different slopes of the surface ten- 
sion coefficient below and above its nil line T = (/i) 
are not unusual since this property is successfully used 
in such well known models as the Fisher droplet model 
(FDM) [33] and the SMM O [37j, but, additionally, 
in the present model the reduced surface tension co- 
efficient is negative (positive) for T above (below) the 
line T — T^(p). As it is argued in [15l [19] there is 
nothing wrong or unphysical with the negative values 
of surface tension coefficient, since in the grand canoni- 
cal ensemble the quantity TE is the surface free en- 
ergy f S urf = e s 

ur f - Ts s ur f of the bag of mean volume 
v, were e sur t and s sur f are the surface energy and en- 
tropy. Therefore, E < means that the surface entropy 
contribution simply exceeds the surface energy part, i.e. 
Ts sur f > e sur f and then f sur f < 0. It can be shown on 
the basis of exactly solvable model of surface deforma- 
tions [35] that negative values of the surface free energy 
is a consequence of very large number of non-spherical 
configurations at high temperatures. To our best knowl- 
edge, the exactly solvable models of the liquid-gas PT 
with negative values of the surface tension coefficient 
provide us with the only physical reason preventing the 
condensation of small droplets into a liquid phase (an in- 
finite droplet) at supercritical temperatures, and, thus, 
they naturally explain the existence of a cross-over both 
in QGP [TH [19] and, probably, in the ordinary liquids 
|39j . For the field-theoretical arguments in favor of the 
negative surface tension of quark gluon bags see [3D[ [JT] . 
Another strong line of arguments in favor of the negative 
surface tension of quark gluon bags at high temperatures 
is provided by the recent analysis of the relation between 
the confining color string tension and the surface tension 
of QGP bag It clearly demonstrates that at the 

cross-over region the surface tension coefficient of large 
bags is unavoidably negative and, as shown in |42j . this 
should lead to an appearance of surfaces with the fractal 
dimension. Note that the fractal surfaces are well known 
in the lattice QCD formulation [43, 44] . but their princi- 
pal role in the lattice entropy maximum formation of the 
confining tube (the so called 'mysterious maximum' |45| ) 
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and their relation to the negative surface tension values 
of such tubes were revealed only recently [42]. 

The phases of the QGBSTM2 include the hadronic 
phase and the QGP which in the fJ> — T plane are sep- 
arated by the nil line of the surface tension coefficient 
and they can be distinguished by the sign of the surface 
tension coefficient (see Fig. 1). For a given \i and tem- 
perature T above (below) the line Ts(/z) the pressure is 
marked by the superscript +(— ) and is defined by the 
equation 



P 



: (T,/i) = T[F H (z ± ,T^) + u(T^)I T (A ± z,^ ± )],(10) 



where the following notations 



P ± (T, M ) 



zm(T,h) are used. The expression (10) for 
+ 



pressures p^ and p~ is determined, respectively by the 
simple poles s + and s~ of the isobaric partition 

The mixed phase of these two phases corresponds to 
the vanishing value of the surface tension and in this 
respect it is similar to the CEP in ordinary liquids [5], 
but, in contrast to the ordinary liquids, the PT in the 
QGBSTM2 is not of the 2-nd order, but of the 1-st or- 
der everywhere at the phase diagram except for the CEP 
where it is, indeed, of the 2-nd order [15]. The corre- 
sponding pressure is given by the essential singularity of 
the isobaric partition ^ 



p M {T,n) = Tz M (T,n). 



(11) 



This equation also demonstrates the meaning of the func- 
tion Tzm{T, fi) which would give the QGP pressure in the 
absence of the surface tension. 

As usual, the pressure of the stable phase is defined 
by the rightmost singularity of the isobar partition ([5| 
(I5HI5] . By construction the PT occurs at T = T 2 (/j) at 
which the simple pole singularities z ± (T, /z) coincide with 
the essential singularity zm(T, fi). The colliding singular- 
ities automatically provide the fulfillment of the Gibbs 
criterion of phase equilibrium. The necessary condition 
for the PT existence is the following relation between the 
parameters of the model 

z M (T, f i) = F H (z M (T,fi),T,ii)+u(T,fi)I T (0,0),(12) 

at T = T c (fi) = Te(jjX 

Assuming that (12) is fulfilled we parameterize the 
shape of the phase coexistence curve T c (fi) in the vicinity 
of the CEP by the constant £ T > 0: 



T, 



cep 



T c{lA ~ (A* - Peep) 



(13) 



The crucial importance of this new index for the QGB- 
STM1 was demonstrated recently [18]. Surprisingly this 
index was not considered in such well known models as 
the FDM [33J and the SMM [H[37], although it is known 
that a related quantity K c introduced in 46 plays a de- 
cisive role in the classification of the CEP stability types 
|47j . Below it will be shown that just this index deter- 
mines the values of the exponent a' . 
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FIG. 1: [Color online] A schematic phase diagram in the plane 
of baryonic chemical potential fi and temperature T. The 
deconfinement PT line T c (/i) is shown by the full curve for 
fi > fj, cep , whereas the cross-over (shown by dashed curve) 
takes place along the line Te(ju) for /i < fi C ep- For fi > fi cep the 
reduced surface tension coefficient E changes the sign exactly 
at the PT line, i.e. T c (ji) = Tgfju) for fi > fi cep . The PT 
region ends at the CEP (filled circle). 



Using the standard definitions for the entropy den- 
sity s and the baryonic density p as T and /i partial 
derivatives of the corresponding pressure one can ex- 
plicitly write the Clapeyron-Clausius equatio n f or pure 
phases ^ = - s p tz s p - \ T=T witn th° h e h> of (10 1. How- 
ever, in the present model there is an additional relation 
for the pressure of the mixed phase pm (ID, w hich by 
construction matches the pressures p~ and p + ( 10 1 for 
the same value of the PT temperature. Therefore, one 



can establish two additional relations of the Clapeyron- 
Clausius type between the partial derivatives of the func- 
tion pm (T, /i) and the partial derivatives of the pressure 
of each pure phase 



dfi c 

It 



sm 
Pm 



t=t c 



T=T C 



where the following notation 

dF H du T 3zm ( dF H 



A, 



dx dx 



It 



dx \ dz 



1 



(14) 



(15) 



for x € {T, fi} is used. 

For further evaluation it is convenient to parameterize 
the behavior of the numerator and denominator in ( 14 ) 



at the CEP vicinity in the same way as it was suggested 
recently in [15] : 



A T 
A 



T=T r 



M I T=T C 



{T cep - T c {p)) x+ e 

(T cep - T c (/i))*, 



(16) 
(17) 



where the finite values of the integral I T (0,0) and the 
functions Fh(z (T,fj,),T,[j,), u(T,fi), Zu(T,fi) along 
with their first derivatives for any finite values of T and fi 
provide the validity of the condition \ > max(0, 1 — ^r). 
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Since the index \ unavoidably appears from an inspec- 
tion of the Clapeyron-Clausius equation, which is a di- 
rect consequence of the Gibbs criterion, then it is quite 
general |18j . Moreover, the parameter \ played an im- 
portant role in separating the different sets of solutions 
for the QGBSTM1 jl8j, but as it will be shown below, 
although such an index appears in the intermediate ex- 
pressions for the analyzed singularities, surprisingly, it 
does not enter any equation for the critical exponents of 
the present model. 



III. THE STANDARD CRITICAL INDICES OF 
THE QGBSTM2 

As usual, the standard set of the critical exponents 
a', /3 and 7' [Sj [13] describes the T-dependence of the 
system near the CEP: 

C p ~ \t\- a \ for t<0 and p = p cep , (18) 



Ap ~ \tf, for t < 0, 
AK T ~ |t|" 7 ', for t < 0, 



(19) 
(20) 



where Ap = (p + — p )t=t c defines the order param- 
eter, C p = ^-(f^)p denotes the specific heat capacity 

at the critical density and AKt = (K^ — K^)t=t c 
is the discontinuity in the isothermal compressibility 
Kt = o(f^)r across the PT line, the variable t is the 

The critical isotherm 



T— T ce 



reduced temperature t 
shape is given by the index 5 [51 [5] 



Pcep - P ~ (Peep - P) 5 for t = 0. 



(21) 



Hereafter the tilde indicates that T = T cep . Note that 
for the phase diagram shown in Fig. [l] one has p = p + 
and p = p + . 

The critical exponent a' describes the T-behavior of 
the specific heat capacity along the critical isochore p = 
p cep inside the mixture of the QGP and hadron phase. 
Similarly to [TH] here we do not use the Yang- Yang for- 
mula [48 to calculate C p , since it leads to a less conve- 
nient representation. Therefore, some details of the heat 
capacity evaluation are given below. We have to note 
that there exist a few suggestions on how to possibly de- 
termine the specific heat capacity at the vicinity of the 
QCD CEP [351 HO] from the experimental data obtained 
in the heavy ion collisions and, hence, the question of 
correct definition of the a' index becomes of crucial im- 
portance for us. 

As usual, the entropy density and the baryonic density 
of the mixed phase are defined via the corresponding val- 
ues of pure phases and the volume fraction of hadronic 
phase A (the volume fraction of the QGP is, respectively, 
1-A): 

p\t=t c = A j0 -| T=Tc + (l-A)p + | T =T c , (22) 
s\ t =t c = A.s-| T=Tc + (1 - A)s+| T=Tc . (23) 



Varying A from to 1 one can describe all states inside the 
mixed phase. Fixing p\t=t c to p cep in (22 1 one can first 
calculate the total T-derivative of the volume fraction A 
along the critical isochore p — p cep and then one can 
use ( 23 ) to determine the specific heat capacity at fixed 



baryonic density 



T 

Pcep 

+ A 



ds+ 
~dT 



dp c dp + 
~dT~dT 

dp c 



dT 



(s- - s + ) + 



dT dT 



(P--P + ) 



T=T a 



T 



ds+ 

It 



dpcdp^ . + _ , d 2 p c 
dT dT + [P Pcep) dT 2 



,(?4) 



T=T C 



where in the second step the Clapeyron-Clausius equa- 



tion ( 14 ) for hadron phase and that one for the QGP 
together with the definition (22) for A were used. Using 



the Clapeyron-Clausius equation ( 14 ) for the QGP one 
can again rewrite ( 24 ) as 



T 

Pcep 

dp. 



ds M 
dT 

d 

dT dT* 
T 



dp c dp M 
dT dT 



+ :^( s+ - s m) 



dT 



Pcep 



(p + 

(Pm 



- Pm) + {p + 

Pcep) 



Pcep) 



ds 



M 



d 2 p c 

(Jj-I- rj-t rj-\ 

dpc dp M ' 



a p c 

dT 2 ^ dT ^ dT dT 



.(25) 



T=T C 



The obtained expression ( 25 ) for the specific heat capac- 



ity coincides with that one obtained for the QGBSTM1 
[18] and, therefore, the a'-index for the models with the 
CEP and with the triCEP are identical. 

In this section we consider the divergent heat capacity 
at the CEP, i.e. a' > 0, while the negative values of this 
index are analyzed in the subsequent section. Then the 
critical exponent a' > describes the temperature be- 
havior of the most singular term in Eq. ( 25 ) . Expanding 



Pm\t=t c into the series of t-powers and using the pa- 
rameterization ( 13 ) of the PT curve one can show that 



the first term staying on the right hand side of ( 25 ) after 
the second equality behaves as £ min ( 1, { T ^ + «' r 2 . Since 
the entropy density and the baryonic density of QGP 
are, respectively, the T and p derivatives of the function 
Tzm (T, p) , which is a regular function of its parameters 
together with its first and second derivatives, then the 



singularity in the second and third terms of ( 25 1 may ap- 



pear from the square of the derivative -^f and it has the 



form (^fr^j ~ i ?T 2 ■ Accounting for these facts, one 
gets the critical exponent a' as 



2 - 2min 1, 



1 



(26) 



This equation shows that a' > for £ T > 1 only, oth- 
erwise a' — 0. As it was mentioned above, the critical 



exponent a' of the QGBSTM2 ( 26 ) exactly coincides with 
that one obtained for the QGBSTM1 18J despite the fact 
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that the phase diagrams of these models have essential 
differences. Moreover, these models have entirely differ- 
ent ranges of the Fisher exponent t: r > 2 within the 
QGBSTM2 while 1< r < 2 within the QGBSTM1. 

To calculate the critical exponents (3 and 7' we need to 
know the behavior of the integrals I T _ 9 (0,0). As it was 
shown in |19) the necessary condition for the 1-st order 
deconfinement PT existence is that such integral should 
be finite for q = 1. This condition provides the range 
of values r > 2 for the Fisher exponent. The validity of 
this necessary condition can be demonstrated as follows. 
Indeed, using the definition of the baryonic density and 
Eqs. ( 10 ) for the phases on two sides of the PT curve 



one can find the baryonic density discontinuity across 
the deconfinement PT line H51: 



Ap 



-T ■ ■ 



(f-f)«U(P) 



^ +uI T ^ (0,0) 



d 

dZ+ \ 
dp ) T _ 



T=T a 



(27) 



T=T C 



for t > 2 since in this case the integrals I T ^^(0,0) and 
7 T _ 1 (0, 0) remain finite. Therefore, the temperature be- 
havior of Ap is defined by the difference of the reduced 
surface tension coefficient partial derivatives calculated 
below and above the PT line. We would like to pay an 
attention to an analysis of the condition Ap > 0. Clearly, 
this condition simply means that the baryonic density 
increases during the deconfinement PT from hadrons to 
QGP, since the baryonic density of latter is higher. 

Since the integral l T - x (0, 0) and the denominator stay- 
ing in the expression after the second equality sign in ( 27 ) 



is positive, then the condition Ap > is provided by the 



inequality 



d^ 
dp 



< 



dp 



which is the case for 



1 



only. Therefore, according to the parameterizations ^ 
and ( 13 1 one obtains 



min(/3 + ,/r), 



for 
for 



>C + 



1, 
1 , 



(28) 



where the notation P ± = + — Jr is used. Since 
in accord with the adopted assumptions the densities of 
pure phases are finite at the CEP, it follows that the 
index (3 is non-negative. 

To find the critical exponent 7' one has to calculate the 
isothermal compressibility Kt = ^(§^)t f° r both pure 
phases. Using the baryonic density definition one can 
rewrite the isothermal compressibility as Kt — ^rf^r- 



Therefore, using Eq. ( 10 ) at the CEP and keeping the 



most singular terms one finds 



a 2 ir 

dp 2 


a 2 s+ 

dp 2 


2 




P 


P+ 2 


d 2 X+ 




dp 2 


dp 2 


P + 


_2 

P 



lW r _ x (0,0) 



1-t£ + u/t-i(0,0) 



T=T a 



(29) 



Such a quantity, as suggested in [13] , is rather convenient 
for the index 7' evaluation. The detailed derivation of the 
expression (29) can be found in the Appendix A. Indeed, 
differentiating Eq. (27 1 with the help of the parameter- 



ization ( 13 ) and recalling that its left hand side is just 
\t\P , one immediately finds 



7 



1 



(30) 



It is evident, that the direct way to calculate this index 
via the isothermal compressibilities Kt of pure phases 
would give us the same result for the present model. 
Thus, at the first glance it seems that the definition (20) 



suggested in [13] is rather appropriate for the case of the 
asymmetric phase diagram, when the 7' indices of two 
pure phases are different. We, however, will return to 
this question again in the subsequent section. 

To calculate the critical exponent S one can use the 
definition of Az at the critical isotherm and get 



P - Peep = T cep (Az + Z M - Z M \cep) 

dz M 



T cep Az + Ap 



dp 



(31) 



cep 1 



where in the second step one has to expand zm in powers 
of Ap = p — p cep and keep the linear term. Similarly one 
can determine the deviation of the baryonic density at 
the critical isotherm from that one at the CEP. Note that 
this critical isotherm is lying outside the mixed phase and 
necessarily it belongs to the high density phase. Then one 
finds 



P Peep 



= T, 



cep 




cep j 



(32) 



Analysis of the functions Az and near the CEP 
shows that it is convenient to substitute the expansion 
/ T (A*s,E) w J T -M-Az/ T _ 1 (^,f)-E/ T _ x (%sf) 
[T5] into Eq. ( |To| for the cross-over states in order to 
obtain 



ApA^| ce p 

1 _ 9F H 

dp 



- Su I T - 



( Az El 
t V 2 ' 2 / 



5 ) 
' 2 I 



(33) 



where it is sufficient to keep only the first order terms 
of the expansion and to use the fact that at the CEP 
Az = 0. According to Eq. (31) the higher order terms 



of expansion could be neglected. 

Using the same steps of derivation one can show that 



T=T C 



8Az _ Ap 
dp 



Mel 



dp 



ul 



1 - 



§Fh 
dp 



2 

s > 
2 



(34) 
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The reduced surface tension coefficient partial deriva- 
tive vanishes at the CEP, whereas its second partial 

derivative diverges at this point. Therefore, from 
Eqs. |33|, |32| and ^ onejinds that_Az - A/i 



and 



as 



g ■ - g . Then using Eqs. (|9|) and (j 13|) one con- 
cludes that 
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(35) 



(301 



Note that the condition S > 1 requires > fl + which 
according to the result found for the index 7' 
consistent with the constrain 7' > 0. 

Despite some similarities the critical exponents of the 
QGBSTM2 are very different from the critical indices of 
the comparable cluster models such as the FDM [33J , the 
SMM 17 and the QGBSTM1 Q3]. Thus, surprisingly, 
the critical exponents of the present model do not depend 
on the Fisher topological exponent r and on the param- 
eter x which relates the mean surface of the QGP bag to 
its volume. It seems to be a unique feature of the QG- 
BSTM2, since the critical exponents of the FDM, SMM 
and QGBSTM1 are expressed in terms these parameters. 
Furthermore, in contrast to the QGBSTM1 critical ex- 
ponents [18], the parameter \ does not appear in any 
final expression for the critical indices of the QGBSTM2. 
This is another surprising fact since just this parameter 
switches the different regimes between the QGBSTM1 
critical exponents [18] . 

Having the explicit expressions for the QGBSTM2 crit- 
ical exponents let us now discuss the question whether 
the obtained results are able to reproduce the indices of 
the 2-dimensional Ising model [51] . of the 3-dimensional 
Ising model [52], of the simple liquids [51] and of the 
0(4) symmetric 3-dimensional spin model [53J [53] (see 
the Table [i I . It is an important question since, as it is 
expected, the universality class of the 3-D Ising model 
coincides with that one of the 3-flavor degenerated QCD, 
whereas the 2+1 flavor QCD falls into universality class 
of the 0(4) spin model [H [55J OH [55J . 

As one can see from the Table [I] 7' > 1 for all discussed 
systems. On the other hand, from the explicit expressions 
of the QGBSTM2 critical exponents one can see that 
within the present model 7' > 1 for £ T < 1 only, which 
immediately leads to a conclusion that a' = 0. Therefore, 
the QGBSTM2 with the traditional definition of index 
a' is able to reproduce the critical exponents of the 2- 
dimensional Ising model only. Such a problem is not a 
new one. For example, the critical exponents of the SMM 
[17j do not reproduce that ones of the simple liquids and 
of the 3-dimensional Ising model since oe' SMM = 0. It is, 
however, believed that such a problem is related to the 
traditional definition of the critical exponent a' . In order 
to elucidate this fact let us study in detail the scaling 
relations for the indices of the present model. 



TABLE I: The critical indices of the 2-dimensional Ising 
model 51, of the simple liquids [5T], of the 3-dimensional 
Ising model [52] and of the 0(4) symmetric spin model 



IV. THE SCALING RELATIONS OF THE 
QGBSTM2 

Now we return to the discussion of the scaling inequal- 
ities ([!]), ([2]) and ([3| of the present model. This inequal- 
ities are well established analytically [5_H7] and experi- 
mentally [SI (5TJ [5J5] . As one can see from the Table II 
these inequalities are exactly established only for an ex- 
act solution of the 2-dimensional Ising model while for 
other systems shown there such inequalities are estab- 
lished within the error bars, which in some cases are not 
very small. From time to time in the literature there ap- 
pear the models [T3J \T7\ [18] and even the experimental 
works (we mean the text-book example [57j discussed in 
the quoted edition of [H]) in which the problems related 
to the Fisher ([I]) and the Griffiths ^ inequalities are 
reported. It is remarkable that the reported problems 
are always related to the inequalities in which the expo- 
nent a' is involved. The situation is somewhat mysteri- 
ous since the formal conditions of the well known Fisher 
theorem [5] proving the validity of these inequalities for 
liquids are fulfilled by the models [T3J [TT] [TS] , but there 
is the range of parameters for which either one of the 
relations ([!]) and ^ or both of them are not obeyed. 
Therefore, it is interesting to verify the scaling relations 
for the indices of the QGBSTM2 which, as we discussed 
in the preceding sections, demonstrate the markable dif- 
ference with that ones of the models [T3J [TT] [TS] . 

The explicit expressions for the critical exponents 
found above allow us to directly examine the scaling in- 
equalities. Despite the usual expectations, the Fisher 
and Griffiths inequalities are not always obeyed, whereas 
the Liberman inequality is fulfilled for any values of the 
model parameters. Let's first demonstrate the validity of 
the Liberman inequality ([3]). Using the explicit expres- 
sions for the indices /?, 7' and S, i.e. Eqs. (28), (30) and 
(35 1, one obtains 



y + /3 (i-£) = i(i-^) 

forr>C + = l 
for C = (+ = 1 ' (36) 



1 - 



min(/3+,/3~) 



/5H 



> 0, 



where an obvious inequality min(/3 + ,/3 ) < f3 + is ac- 
counted for. From the Liberman scaling law ^ one im- 
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a' + j3{S+l) 
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2.000412(5000) 


1.9425(55) 


2.02(9) 


i+P(l-5) 





-0.000052(2000) 


0.0775(212) 


-0.01(6) 



with all its derivatives up to the second one are finite, 
we c onclude that the first term on the right hand side of 

2)+/3 



(1391) behaves as £ min(0 'e? 



Similarly one finds that 



the second term on the right hand side of (39) behaves 

_J 2+8 

as t ( T . This analysis shows that the first of the two 
discussed terms is the leading one and, hence, we have 



TABLE II: Scaling relations between the critical exponents 
taken from the Table[l] The uncertainties were calculated from 
their values given in the Table [T] using the error determination 
method for indirect measurements 1581. 



mediately deduces the relation between the scaling in- 
equalities @ and @ a' + fi(6 + 1) < a' + 2fi + 7'. Then 
using the explicit expressions for the QGBSTM2 critical 
exponents one can get the following result for the Fisher 
and Griffiths inequalities 



a' + 13(6 + 1) < a' + 2(3 + 7' 
= -2 (min(l, i) - 1 - - 
= 2 max(l, Jr) — 7' . 



1 



(37) 



Note that in the evaluation of the second equality above 
we employed the expression for index 7' and the following 

0-1- rr = min(-l,-7^) = 



sequence of steps min(l, 
— max(l,|! r ). Eq. (p7k gives the estimate from above 
for the left hand side of the Griffiths inequality ^ and 
also it demonstrates that the scaling laws ([!]) and ^ 
may be broken for the QGBSTM2, if 2 max(l, p) - i < 
2. Moreover, Eq. (37) clearly shows that, if the Fisher 



scaling inequality is fulfilled, this does not guaranty that 
the Griffiths one is obeyed. One of the consequences 
of Eq. (37) is an appearance of the following relation 



between the QGBSTM2 critical exponents 
a' + 2(3 + 2 7 ' > 2 , 



(38) 



which is based on an obvious inequality 2 max (l, ...) > 2. 
Note also that for £ T > 1 the inequality ( f38| turns into 
an equality since max(l, ^V) = 1 in this case. 

In order to recover the scaling inequalities |l]) and ^ 
we follow the suggestion of [13] and replace the index 
a' by a' s , where the latter describes the temperature de- 
pendence of the specific heat capacity difference of two 
phases at the CEP, i.e. AC = (C p + - C p - ) t =t c ■ The 
idea behind such a suggestion is to get the most singu- 
lar term from the difference of the specific heats capacity 
of two pure phases. Then from the Clapeyron-Clausius 
equation one can find that 



AC = T c - 



p+p- 



dT 
p" 



sm 



dfi c 
~dT 



Pm 



rr.d^C d 

Tc lTdT ln p 



(39) 



~ t" in the vicinity of the CEP, then using 
(13) and the fact that the function zm(T,p) together 



Since p" 1 



1 

w 



1 



(40) 



Note that a' s > for A < 2 - j3 only. Using the a' s 



expression ( 40 1 one can get 



a' s + 2/3- 



7 



>2, 



(41) 



which is in a complete agreement with the expectation of 
|13) . However, such a replacement does not recover the 
Griffiths inequality. Indeed, in this case one finds 



■13(1 + 6) 



1 



(3 



(42) 



Obviously, the right hand side of ( 42 1 is less than 2 for 
/9 = P~ < /3+min(l,2£ T ). Thus, the QGBSTM2 gives 
the first and explicit example that the definition of the 
a' s index does not recover the Griffiths scaling inequality, 
although it, indeed, redeems the Fisher scaling law. The 
reason that the a' s definition worked well in all preceding 
models, but fails for this one, is that within the mod- 
els H3 HI] the PT corresponds to the change of the 
leading singularity from the simple pole of the isobaric 
partition partition to its essential singularity, whereas in 
the QGBSTM2 the leading singularity above and below 
PT is the simple pole. Then for the models [HEIDI] the 
divergent terms in the difference AC = (C p + — C p - )t=t c 
did not cancel each other because they have rather dif- 
ferent analytical behavior, while in the QGBSTM2 the 
most divergent terms in the specific heat capacity of two 
phases coincide with each other and they simply can- 
cel each other in the expression for AC. Therefore, the 
QGBSTM2 provides us with a direct evidence that such 
a simple and attractive definition of the a' s index [13] 
which, so far, was designed for the asymmetric phase 
diagrams and worked well for the models [131 El ITS] . 
should be replaced now by a proper one. 

Keeping in mind these facts we introduce a new def- 
inition for the index a' which is based on the behavior 
of the linear combination of specific heat capacity of two 
phases C tot = [K+C+ + K~ C-] T=Tc ~ \t\~< with the 
positive and nonvanishing coefficients K ± > which in 
general may depend on T or [i. The inequality > 
guarantees that no term in C tot is missing or cancelled. 
On the other hand the condition K ± \t=t c > provides 
that the index a' c has the maximal value among that one 
of pure phases. Note that this condition leads to a differ- 
ence with the traditional definition of the index a' since 
the coefficients K + \t=t and K~\t=t are not related to 
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the fractions of pure phases. Such a property allows us 
to avoid the situation, when the curve at which the spe- 
cific heat capacity is found does not match the boundary 
with one of the phases which was the case for the models 
analyzed in [T7] . 

The index a' c definition has a clear physical interpre- 
tation: when approaching the CEP the density fluctu- 
ations get so strong that the specific heat capacity of 
both phases contribute into the measurable value and, 
hence, the largest term determines the critical exponent 
a' c . Then the coefficients K + \t=t c and K~\t = t c define 
the weight of corresponding pure phase into the specific 
heat capacity of their mixture. Obviously, the maximiza- 
tion of the a'-value is not well suited for the pure phases 
with asymmetric properties, but also it should increase 
the magnitude of the left hand side of the Fisher and 
Griffiths inequalities and, thus, it should weaken the nec- 
essary conditions to obey them. 

Let's find out the index a' c . Making the same steps as 
in deriving Eq. ( 39 ) we obtain 



Ctot — T 

- T 
K 



K~ 
P" 



dS M : 

dT It = t - 
dfj, c 



K+ d ( + 
-p^dT\ {p - pM) dT 



d 

dT 



(P - Pm )) 



dpc 
~dT 



(43) 



T=T C 



From this formula one can draw an important conclusion 
that negative values of the index a' c are possible, if and 
only if each full T-derivative on the right hand side of 
(|43| including ^1 and _ vanish at the CEP. 

° al \cep al \l—l c 

The necessary conditions for a' c < are -£f-\___= and 



ds 



M 1 



dT 



lT=T c 



i*r> 



(44) 



(441 



for a' c > 0. The 
is that the function 



where uj > for a' c < and uj 
geometrical meaning of Eq 
sm{T, Hcep) has a kink point at the critical temperature 
for uj > 0. The case uj = does not add anything new 
to the above evaluation of the indices a' and a' s . Us- 
ing the parameter uj one ca n sh ow that the first term on 



(|40b behaves as t^^'e 



the right hand side of Eq. 
Since the coefficients K^\t=t c are positive, then each 

term Up* - Pm)^H in (44) snomd vanish in- 

dependently. As shown in the Appendix B, employing 
the parameterizations ([9| and (17) one can find that 



max[(p + — pm), (p~ — Pm)} ~ t p at the vicinity of the 
CEP. Taking into account this result one can deduce that 
the term staying i nsid e the square brackets on the right 

I I J_ + fl_2 

hand side of Eq. (|43|) behaves as t« T 



Thus, this 



analysis allows us to determine the critical exponent a' c 
as 



1 



+ P-uj,2) -p- 



1 



(45) 



Similarly the exponents a' and a' s can be found for the 
nonzero values of to 



a'lujto = 2 - 2min 1 



1 U 1+(T 

— l-i L_ 

e T ' 2' 2 



'i.L-n = max | - u},2\ - /3 - ~ 



, (46) 
(47) 



Comparing Eqs. (45) and (47) one can easily see that 
a' > a' for the same value of uj. Moreover, for the case 



uj > 2 one finds that a' = a' 



the 



index a' c can be essentially larger then a' s since ft > 0. 
Therefore, it is evident that in the models [131 [HI HI] 
in which the a' s index recovers the scaling inequalities 
([l])-(|3]) the a' c index should recover them too, since by 
the construction it is not smaller than a' sl but for other 
models the properties of this index should be studied. 

Let us demonstrate explicitly that the proposed defi- 
nition of the critical exponent a' c is more adequate with 
respect to the QGBSTM2 scaling inequalities. From the 
expression (30) for the index 7' one gets ^ = 7' + p. 



Substituting this result into ( 45 ) one finds 



-uj , for uj < 7' + 2/3 -2 

2 - 2(3 - i , for uj > i + 2(3 - 2 = 



(48) 



and, hence, a' c > 2 — 2/3 — 7' which is just the Fisher 
inequality written for the oi c exponent. Note that in con- 
trast to the a' c definition the Fisher inequality for the 
a' s definition is not fulfilled anymore. For instance, if 
a' c = —uj, then from Eqs. (47) and (|48| it follows that for 



the same value of the index uj the relation a' s +3/3+7' > 2 
holds which, obviously, differs from the Fisher inequal- 
ity. The considered example clearly demonstrates that 
for a' < the most general definition of the index a' s 
leads to the difficulties with the Fisher scaling relation 
and, hence, it has to be replaced by the a' definition. 



Since the Liberman inequality ( 36 ) does not depend 
on the a' definition, it holds, but here we would like 
to pay an attention to the case when it turns into an 
equality. From Eq. (36) one concludes that this occurs 
for (3 = (3 + , which applied to the relations for the indices 
6 and 7 yields a very important relation 



y + P = (36= -L 



(49) 



that should hold for any model in which the scaling re- 
lation ([3 1 is obeyed as equality, for instance, for 0(2), 
0(3) and 0(4) spin models. For such models the expo- 
nent £ T = ^ is the ratio of the so-called magnetic yh 
and thermal y t exponents [SHI [Ml ISO] and it defines the 
curvature radius of the phase diagram in ji — T plane 
|59j . If, however, the scaling law ([3| is obeyed as an in- 
equality, the relation (49) reads as an inequality too, i.e. 

(3se<i- 
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For the Griffiths inequality J2| the direct substitution 
of the expressions (28), (30), (35) and (45) yields 



+ p(i + s) = 



+ p(l + 6), fbrw</? + 
i(l-^),fora;>/? + 



2, 
2. 



(50) 



Note that the right hand side of Eq. ( 50 ) can be smaller 
than 2 in two cases. Indeed, if the index oj obeys the 
inequalities /3(1 + 5) — 2<oj<f3+^ T — 2, then the 
Griffiths scaling law (2| is broken down. The other pos- 
sibility for the break down of this scaling law corresponds 
to the lower equality in ([501, i.e. for 



(51) 



a' c + p(l + 6) = 2 - 5(/3+ - p) , 



that occurs for a; > (3+^—2 and P < P + . Thus, 
from these examples one might conclude that even the 
improved definition ( 45 1 of the index a' which perfectly 



works for the Fisher inequality does not save the Griffiths 
scaling law. We, however, would like to stress that these 
examples clearly show us that the problem of formulat- 
ing the scaling inequalities for the asymmetric phase di- 
agrams is much deeper and cannot be resolved by the 
redefinition of the indices. It also requires the modifi- 
cation of traditional scaling relations. Indeed, using the 
Fisher inequality with the a' c exponent and the explicit 
expressions for the indices 7 and 5 we obtain another 
inequality 



a' c + P + p + S>2. 



(52) 



which is an analog of the Griffiths inequality for the QG- 
BSTM2 exponents. 



Let us show that the scaling law (52) cannot be bro- 
ken within the QGBSTM2. For this purpose we assume 



that a' c + P + P+5 = a! c + P + 



< 2. Then, 



paring the latter assumption with the upper equality in 
(50), we arrive at the contradiction, since in this case 
uj = — a' c < P + jr — 2, but our assumption is equivalent 

to the different inequality, i.e. ui 



-ct c > P + i - 2. 



Analogously, from the lower expression in ( 50 1 it follows 
that a' c + P + P + 5 = 2 and we again obtain the contradic- 
tion with the original assumption which, as we proved, is 
the false one. 

From Eq. (50) one can deduce that Eqs. (49) and (51) 



hold as equalities for w > /3 + — 2 and = P + , and, 
hence, the Fisher scaling relation ([I]) becomes an equality 
too. Thus, in this case all scaling laws ([l])~([3]) are fulfilled 
as equalities and this is the traditional scaling regime. In 
all other cases defined by Eq. ( 50 ), i.e. for P(l + 5) — 2 < 



uj < ,5+^-2 or for w > /3+|^-2and^ < the Fisher 
and Liberman scaling laws are obeyed as inequalities, but 
the Griffiths one in its usual form ^ is broken down 



and, hence, we suggest to use the inequality ( 52 ) instead 
of ([2]). Indeed, the latter, as we proved above, is valid 
in more general case than the original Griffiths one and 



also it seems to be more natural, since, in contrast to pj), 
for P + > P = p~ the high density phase exponent d in 



( 52 ) is multiplied by other high density phase expo nent 
P + . Since in this regime the generalized inequality ( 52 ) 
is valid, we call it the generalized scaling regime. 

Also we have to stress that the break down of the tradi- 
tional Griffiths inequality ([2]) does not require some spe- 
cial conditions. Indeed, it is easy to see that, if the Fisher 
relation ([!]) is fulfilled as equality and the Liberman one 
([3]) is obeyed as a strict inequality, then the traditional 
Griffiths inequality ^ would be broken, which is not 
the case for the inequality (52). Note that for the QG- 



BSTM2 exponents the Liberman inequality can be also 
generalized similarly to the Griffiths inequality (52). In- 



deed, rewriting the first equality in the expression ( 36 ) 
one finds 



7' + P - p + S = . 



(53) 



which is stronger than the usual Liberman scaling in- 
equality, since it is possible that P + > p. Here it is 
appropriate to discuss the definition ( 20 ) of the critical 
index 7' suggested in [T3]. Similarly to the definition of 
the a' s index [13] such a definition is simple and it may 
provide one with the maximal value of the compressibil- 
ity exponent of two pure phases. However, it has the very 
same defect as the a' s definition, namely it does not guar- 
anty that for some models with the asymmetric phase 
diagrams the leading terms of two compressibilities may 
simply cancel each other as we found such a possibility 
for the QGBSTM2 a' s index and this would lead to the 
problems with the scaling relations. Therefore, in order 
to avoid these problems for such cases we suggest to em- 
ploy the definitio n of the 7' exponent which is similar to 



the a' c definition ( 43 ) in which the heat capacities C are 



replaced by the corresponding thermal compressibilities 



k + K+ 



t 



(54) 



with the positive and nonvanishing coefficients > 0. 
The interpretation of such coefficients is similar to that 
one introduced for the a' c index. 

After the thorough discussion of the QGBSTM2 criti- 
cal exponents it is not surprising that for the traditional 
scaling regime the present model with the a' c and ui def- 
initions is able to reproduce the critical indices of the 
2-dimcnsional Ising model, of the simple liquids, of the 3- 
dimensional Ising model and of the 0(4) symmetric spin 
model for which, as one can see from the Table II, the 
inequalities ([l])-([3]) are well established. The list of the 
corresponding QGBSTM2 parameters is given in the Ta- 
ble III. However, an existence of the generalized scaling 
regime for the QGBSTM2 with the most favorable defi- 
nition of the a' c index for which the Griffiths inequality 
is broken down is rather surprising, since it is widely ac- 
cepted that the validity of the Fisher and the Griffiths 
([2]) relations does not require any additional assumptions 
except for the conditions of the Fisher theorem J5[ [H [SB] . 
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2D Ising 


3D Ising 


Simple 


0(4) 




Model 


Model 


Liquids 


Model 


s 


8 
15 


0.6010(1) 


0.631(2) 


0.55(1) 


c+ 


1 


1 


1 


1 


c 


1 |> 1 


1 


1 


1 > 1 


e 


1 


1.0112(2) 


0.937(23) 


1.20(4) 


r 


1 |>0 


0.9903(2) 


0.910(25) 


1.20(4)| >0 


UJ* 











0.19(6) 


UJ** 


> o 


> o 


> o 


> o 


UJ*** 


> o 


> o 


> o 


> 0.19(6) 



TABLE III: The QGBSTM2 parameters which describe the 
corresponding exponents given in the Table [j] The values of 
the parameter uj extracted from the expressions for a' , a' g and 
a' c are marked with *, ** and ***, respectively. The symbol 
means that it is impossible to find the value of the corre- 
sponding parameter which allows us to describe the critical 
exponents shown in the Table [JJ 



Now the QGBSTM2 provides us with an explicit example 
that this may not be the case for the Griffiths inequality 
^ and the latter has to be replaced by (52 1. 

It is worth to note that from the physical point of view 
the successful description of the critical exponents be- 
longing to the different universality classes by a single 
model not only shows that the QGBSTM2 is very general, 
but it also evidences for the fact that the physical clusters 
or bags employed in this model are, indeed, the relevant 
degrees of freedom at the CEP for all analyzed systems, 
despite the very different physics exhibiting by the orig- 
inal Hamiltonians of the corresponding spin models. On 
the one hand this conclusion completes the finding [ST] 
that the FDM correctly describes the distribution of 
large clusters of the 2- and 3- dimensional Ising model in 
the wide range of temperatures, including the CEP. On 
the other hand the QGBSTM2 generalizes the results of 
the noninteracting cluster models which quite 

successfully describe the CEP properties of real liquids. 
Since the QGBSTM1 [TS] is also able to reproduce the 
critical exponents [IB] of the same universality classes as 
the QGBSTM2, we may hope that these models can be 
effectively used to describe the QCD endpoint properties 
and they can help to experimentally distinguish the CEP 
case from the tricritical endpoint case. 



V. CONCLUSIONS 

The practical necessity to describe the thermodynamic 
functions of the QGP with the (tri) critical endpoint 
which has the required properties stimulated the develop- 
ment of a variety of the exactly solvable statistical models 
[T51 [MHSE] . Since it is not exactly know whether 

the QCD phase diagram endpoint is critical or tricritical 
it was necessary to developed the exactly solvable mod- 
els for both cases. However, in contrast to the tricritical 
endpoint case which was worked out in many versions 
[Tfjl UM [T51 IrjlrlSg] . a formulation of a realistic exactly 



solvable model with the CEP required many years and 
took many additional efforts, because neither the phys- 
ical mechanism of its generation nor the mathematical 
properties of such a model were known before its devel- 
opment in [19j . The main result of |19j clearly demon- 
strates that the traditional framework of liquid-gas PT 
models [SJ [TU [33J [37J has to be essentially modified in 
order to degenerate the 1st order deconfinement PT into 
a CEP and into a cross-over. Furthermore, on the one 
hand the liquid-gas PT models require the vanishing of 
surface tension coefficient at the CEP, and on the other 
hand the positive values of the confining tube entropy 
unavoidably demand for the negative surface tension co- 
efficient at the cross-over region [42j 69|. The natural 
conclusion that follows from these findings is that at the 
cross-over region there should exist the curve of nil values 
of the surface tension coefficient which must pass through 
the endpoint of the 1st order deconfinement PT. If the 
nil line of the surface tension coefficient does cross the 
deconfinement PT curve at the endpoint only, then, as 
shown in |15j , such an endpoint of the deconfinement PT 
is the tricritical one and this is the case of QGBSTM1. 
If, however, the nil curve of the surface tension coefficient 
matches the deconfinement PT curve as described above 
(see Fig. 1 for details), then the model has the CEP and 
it corresponds to the QGBSTM2 [15] . 

The present paper is devoted to the analysis of the 
critical exponents of the QGBSTM2 which are necessary 
to be studied before applying the present model to de- 
scribe either the lattice QCD thermodynamics or various 
experimental data. Such an analysis allowed us to find 
some general restrictions on the model parameters and 
figure out the important relations between them. Also it 
is found that the QGBSTM2 exponents essentially dif- 
fer from that ones of the FDM, SMM and QGBSTM1 
despite many similarities between these models. Thus, 
in the FDM and the SMM the critical exponents can be 
expressed in terms of a few major input parameters like 
the Fisher topological exponent r and the exponent k 
which is usually related to a dimension, whereas the QG- 
BSTM1 critical exponents also depend on two recently 
introduced parameters £ T and \ [15] . In addition to all 
these parameters the QGBSTM2 depends on the input 
exponents £ , £ and uj, but in this work it is found that, 
in contrast to the FDM, SMM and QGBSTM1, the QG- 
BSTM2 exponents do not depend on the input exponents 
h, x an d t, although here we for the first time found that 
the present model formulation is valid for r > 3 only. 

Here we also showed that neither with the tradition- 



ally defined index a' of ( 18 ) nor with the specially intro- 
duced one a' s [J3] it is possible to fulfill the Fisher and 
Griffiths scaling inequalities, while the Liberman one is 
always obeyed. Such a situation is known for the mod- 
els [T5rfTrj] . but in the case of present model the index 
a' s does not in general allow one to recover the scaling 
inequalities ([I]) and ^ as well. However, in the present 
work we proved that for a physically motivated definition 
of the index a' c which corresponds to the maximal value 



12 



of the index a' of two pure phases, the QGBSTM2 expo- 
nents recover the Fisher scaling inequality 0, while the 
Griffiths inequality ^ in its traditional form is fulfilled 
only in a case, when all scaling laws ([l])-((3]) become the 
equalities. It is also shown that for this traditional scaling 
regime the QGBSTM2 exponents reproduce the values of 
critical indices of the 2-dimensional Ising model, of the 
simple liquids, of the 3-dimensional Ising model and of 
the 0(4) symmetric spin model and this evidences for the 
fact that the bags (or physical clusters) used in the QG- 
BSTM2 are the relevant degrees of freedom at the CEP 
for all studied systems despite quite a different physics 
of their phases. The present analysis demonstrates that 
the negative values of the a' c index can be achieved with 
the help of newly introduced exponent uj in ( 44 ) , which, 



as we explicitly demonstrated, plays a decisive role in the 
recovering of scaling inequalities. 

Besides the traditional scaling regime, we found that 
the QGBSTM2 exponents have the generalized scaling 
regime, for which the Griffiths inequality should be gen- 
eralized to Eq. (52 1. Also here we showed that for the 



QGBSTM2 critical indices the Liberman inequality can 
be similarly generalized to a stronger relation (53). It is 



clear that a similar results should be obtained for the sys- 
tems whose surface tension at the vicinity of the CEP can 
be reduced to the form of Eq. ^ . The case of the asym- 
metric phase diagrams with other surface tension param- 
eterization should be analyzed both experimentally and 
theoretically. The example of the present model suggests 
that for the asymmetric phase diagrams the Fisher scal- 
ing inequality can be recovered by the proper definition of 
the index a' , while for such phase diagrams the Griffiths 
scaling inequality in its usual form can be broken down. 
This can occur for the linear temperature dependency of 
the surface tension coefficient in both pure phases since 
£ + = £~ = 1 which is not too restrictive as it is demon- 
strated by the FDM and by the exactly solvable model 
of surface deformations 38 . The other important con- 
straint on such systems is that their surface tension coef- 
ficient should change its sign while crossing the PT line. 
Clearly, such a property is rather unusual for the ordi- 
nary liquids and is more resembling the transition from 
a mixture of gases into a strongly interacting plasma, 
which, as one can argue P21 H0H12"] , is the case for the 
QGP. 

Also we believe that for the asymmetric phase dia- 
grams the question of the proper definition of the critical 
exponents and the scaling relations between them should 
be thoroughly investigated both experimentally and the- 
oretically In particular, in many works the index a' is 
difficult to determine directly and, hence, it is usually 
found from the hyperscaling relation [8 , 56 or even from 
one of the scaling equalities or ([2]), while, as it is sug- 
gested by the above analysis, in some cases one can get 
the unexpected results for the traditional scaling inequal- 
ities. 
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VI. APPENDICES 

A. The isothermal compressibility 

In this appendix we give the details of the isothermal 
compressibility evaluation. As usual, one can simplify 
the isothermal compressibility as Kt = ^^| T = ^"f^- 
An explicit expression for the baryonic density of pure 
phases is as follows 



T 



dz M 



A, 



1 - + Ulr-l 

from which one directly obtains its fi derivative as 
~d 2 z M 



dp ± 
dfi 



T 



+ h + h 



(55) 



(56) 



where the following notations are introduced 



dAy, 

dfi 



dfi 2 al T-X Qfl dfJ, 



dfj. u \ dfi 



d£± j 
dfi j t-2* 



(57) 



and 



A, 



dfi 



dF H 



uI T - 



d 2 F 1 



H 



du 



dsdfi dfi 



/dA ± z T as± T 

+ U I — / r _ 2 + -7{ 



V dfi 



dfi 



(58) 



All the regular terms in the expressions (56), (57) and 



(58) can be neglected, since the isothermal compressibil- 



ity behavior near the CEP is defined by the singular ones. 
Hence, keeping only the singular terms we get 



T- 



dp ± 
dfi 

A„ 



dfj, \ dfj, 



dfi 



Ir-2 



1 



dF H 
ds 



uI T -i 



dA ± z t , as ± r 

2 r-2 + -ajrU- 



1- 



ds + U1 T 



-ul r - 



ds 



uI T - 



(59) 



In order to calculate the critical exponent 7' one, ac- 
cording to the suggestion of [T3], has to evaluate the 
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difference of the isothermai compressibilities across the 



PT line AK T 



dp- 
dp 



T=T a 



Requiring that 



1 dp+ 1 

p+ dp p- 

such a quantity remains finite everywhere at the phase 
coexistence curve, except for the CEP, one has to de- 
mand that all the integrals I T - q in Eq. ( 59 ) are finite for 
T = T c , including T = T cep . The analysis shows that it 
is sufficient to require that all such integrals, including 
Jr-2, converge for r > 3 at the PT line, where Az = 
and E ± = 0. Then, the only possible singularity at the 
CEP in the expression for -#^- is provided by the term 

which is proportional to a a s 2 . Then for the isothermal 
compressibility one finds 



dp 2 



±\-2 



2V± 



d 2 E 
dp 2 



TuI T - 



1 



dF B 
ds 



uIt 



(60) 



Just this expression is used to evaluate the right hand 
side of Eq. (pi. 



B. Volume fraction of pure phases 

In order to establish the relation between the indices 
(3 and x h is convenient to analyze the behavior of the 
pure phase volume fractions at the critical isochore near 
the CEP. The QGBSTM2 critical isochore completely lies 
inside the mixed phase. The critical baryonic density can 
be described by the volume fraction of hadrons A (the 
QGP volume fraction is 1 — A, respectively) as 



Peep = Ap |t=t c + (1 - A)/9 + | T=Tc , 



(61) 



where p~ and p + are, respectively, the baryonic densities 
of hadronic phase and QGP. From the previous equation 
one finds 



A 



Peep 



P + - p- 



and 1 — A 



Peep P 



P~ 



(62) 



Using an evident relation p cep = pu\cep and the defini- 
tion of the baryonic density one obtains 

dz M 



(P d 



Pcep)T=T c 

A - 25± 

P' dp 



dp 



cep 

ulr-x 8zm 9zm 



1 



dF„ 

dz 



uI T - 



dp dp 



cep 



(63) 



T=T C 



Applying the parameterization defined by Eqs. Q, (13) 
and (17 1 to the right hand side of (63), one can easily 



t 



min(x,/3 : 



demonstrate that either (p ± — p C ep)T=T c 

for C ± = 1 Or (p± - Peep)T=T c ~ t****^ for C± > 1. 

Note that for this proof it is necessary to expand the 



difference 

op 



dz M 
dp 



cep 



up to linear terms of p—p ce t and 
T c — T cep . Finally, using the definition of the index /3, i.e. 
(p + — p~)t=t c ~ t 13 , and the feature of the present model 
( + = 1 wc find the following behavior of the volume 
fractions 



A 

1-A 



min( x ,/3 + ,l, i ^)-/3 
Lmin(x,l,^)-/? ) for c - 



= 1, 
> 1. 



(64) 
(65) 



in the vicinity of the CEP. Combining these re- 
sults with the fact that both A and (1 — A) are fi- 
nite at the CEP, we obtain a very important con- 
sequences that min(x, f3 + , 1, — j3 > 0, and ei- 
ther min(x, /?~, 1, ^r)(-=i — ft > for (~ = 1 or 
min(x, 1, pr)f->i — (3 > for Q~ > 1. An explicit ex- 



pression for the index j3 ( 28 1 allows us to rewrite all these 
conditions as 



j3 < min I 1 



(66) 



which holds always. From the inequality above one im- 
mediately deduces that f3 < x- These results are im- 
portant for finding the index a' c . Indeed, recalling the 



expression ( 28 1 one can write 



min min \ x, (3+ , 1, , min \ x, /3 , 1, -^j 
for the case (~ = ( + = 1 and 



X,/3 + ,E 



min x, 1 



= P, 
(67) 

/3, (68) 



for the case (, > ( + = 1. From the above equations 
one can immediately deduce that max((p + — Pm), (p~ — 
Pm)) ~ at the vicinity of the CEP, since /? > 0. Thus, 
here we showed that the term staying inside the square 
brackets on the right hand side of Eq. ( 43 ) , indeed, be- 
haves as t « T +l3 2 . 
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